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The general solution of the inverse Frobenius-Perron problem considering the construction of 
a fully chaotic dynamical system with given invariant density is obtained within the class of one- 
dimensional unimodal maps. Some interesting connections between this solution and the approach 
via conjugation transformations are eluminated. The developed method is applied to obtain a 
wide class of maps having as invariant density the two-parametric beta-probability density function. 
Varying the parameters of the density a rich variety of dynamics is observed. Observables like 
autocorrelation functions, power spectra and Lyapunov exponents are calculated for representatives 
of this family of maps and some theoretical predictions concerning the decay of correlations are 
tested. 



I. INTRODUCTION 

One-dimensional iterative maps represent a very useful tool to understand the physics of complex nonlinear systems. 
In particular they can be used to investigate the different routes of nonlinear dynamical systems from regular to chaotic 
behaviour to simulate physical systems showing anomalous diffusion Q or to analyse the underlying dynamics 
in time series with j noise in their power spectrum . In many cases and in particular from an experimental point 
of view the only information available for a physical system are its statistical properties (ergodic measure or time 
correlation function) and one is asked to design a consistent dynamical law. In a recent paper JjJ we examined the 
possibility to construct an one dimensional discrete dynamical system with a given invariant density. This is the so 
called inverse Frobenius-Perron problem (IFPP) |p|-p^[|. 

It could be shown that using certain symmetry requirements the IFPP possesses a unique solution. We applied our 
method to calculate numerically the symmetric unimodal maps with the invariant density given by the symmetric 
beta distribution. However the solution presented in is a partial and very specialized one. Depending on the chosen 
symmetry constraints we can get an infinite number of such solutions. In the present paper we go one step further 
and present the general solution to the IFPP within the class of unimodal maps. A central role is hereby played by 
the function hf which parametrizes the space of the solutions to the IFPP. This solution enables us to explore in 
depth the connection between the dynamics and statistics in one-dimensional systems. Our investigations concentrate 
on the wide family of the beta distributions. The IFPP is solved for the asymmetric case when the two exponents 
in the expression for the density are unequal. From the general solution we obtain two special solutions using two 
different forms for the function hf. We then calculate the time-correlation functions and the power spectra of the 
corresponding dynamical systems. The question of the dependence of the statistical characteristics (singularities of 
the density, long time behaviour of the autocorrelation function) of these systems on their dynamical characteristics 
(critical points) is addressed. Next we compare for a given density the dependence of the statistical properties on 
the form of the function hf. A large variety of the statistical and dynamical behaviour is obtained. Some theoretical 
predictions concerning the asymptotic decay of the correlations or the value of the Lyapunov exponents fl4|-|l"7j| in 
fully chaotic dynamical systems can be tested within the class of our examples. 

The paper is organized as follows: In section II we present the general solution to the IFPP within the class of 
unimodal maps. In section III we derive from the general solution of section II a wide family of unimodal maps 
having as invariant density the beta distribution. This family of maps corresponds to two special solutions of the 
IFPP: the symmetric betamap (SB) and the special asymmetric betamap (SAB). In section IV we discuss the main 
characteristics of the maps belonging to the beta family and we study the rich dynamical behaviour which they 
provide. In section V we calculate the autocorrelation function, the power spectrum and the Liapunov exponent 
for representative members of the beta family. A connection between the dynamics and the asymptotic decay of 
correlations is established. Some theoretical models explaining this behaviour are tested. Finally section VI contains 
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the conclusions and gives a brief outlook with respect to the inclusion of the correlation behaviour into the inverse 
Frobenius-Perron problem. 

II. GENERAL SOLUTION OF THE IFPP FOR UNIMODAL ERGODIC AND CHAOTIC MAPS 

Let us consider the variety of dynamical systems which belong to an arbitrary but fixed invariant density. The 
solution to this inverse problem is of great interest for the numerical simulation of real physical systems as well as the 
understanding of the relationship between the functional form of the map and the statistical features of its resulting 
dynamics. In the present section we establish a general representation of all ergodic and chaotic unimodal maps 
with a given invariant density. This corresponds to the general representation of the solution of the so called inverse 
Frobenius-Perron problem (IFPP) [p[-p^[ within this class of maps. The starting-equation for the construction of the 
map is the Frobenius-Perron equation: 

P(y)\dy\= p{xi)\dxi\ (1) 

where the summation runs over all preimages of y (for unimodal maps i = L(left), R(right)). For any given unimodal 
f(x) the right or left preimage of y is determined, once the corresponding other one is given. The essential feature of 
our approach is the inversion of this expression: A prescribed relation between the two preimages reduces the number 
of independent differentials on the rhs of eq.(|l|) and allows a following integration. Such a relation may be described 
by a function h / (x) , mapping the left preimage onto the right one: 

hf . [0, Xmax] * [^niaxi 1] 

xr = hf(xh) with f(x L ) = /{xr) (2) 

where x max is the position of the maximum of the unimodal map. hf(x) is a function monotonously decreasing on 
the defining interval and it is differentiable with the exception of a finite number of points. It obeys the equations: 

h' f (x) < x G [0, 

hf(0) = 1 hf(x max ) = x max (3) 

Substituting hf in eq.(|l]) yields: 

p(y)dy = p(x L )dx L - p{hf(xL))h' f {xL)dxL (4) 

Using the definition p(x) — pit) dt for the invariant measure eq.(|j) can be integrated in < x < x max providing 
us with the left part /l(x) = /(^)|[o,x ma x] °f the map f(x): 

f L (x)= t i- 1 [(i(x)-n(hf(x))+l] (5) 

with p(x) being invertible and normalized, i.e. p(l) = 1. The right part of the map, fn{x) = f (x)\\ Xmttlt \-\, is obtained 
by substitution of x — > hj 1 (x) in eq.(^J). It follows that the map f(x) parametrised through the function hf(x) and 
characterised by the invariant measure fJ-(x), is given by: 

^i^WV^-^ + i] ; x max <x<i\ = 1 -^-^ H f^ ^ 

where Hf(x) is given by: 

H fW-\hf(x) ; z max <*<l (7) 

With the above eq. (|J) we have found an interesting representation and answer to the IFPP which is the central topic 
of this paper: All unimodal maps with prescribed invariant measure are given by (^|), hf{x) taking on all possible 
functional forms obeying (^|) and (^|). This shows that a fixed invariant measure is a relatively weak constraint to a 
solution of the inverse problem and there is still a considerable freedom to model the mapping. 

Before we apply our method to analyse specific examples we mention an interesting connection between this approach 
and the approach via conjugation transformations which can also be applied to construct a map with an arbitrary 
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given (in general nonsymmetric) invariant density. To our knowledge the conjugation transformation method has 
been applied in the literature only to the restricted case of the construction of a symmetric map with prescribed 



invariant density [[15 16 1. We will show below how this approach can be extented to include all unimodal maps. First 
we recall the basic features of the conjugation transformation approach. More details can be found in the literature 
p8| , ^9| . Using as a basis the tent map given by: t(x) = 1 — \2x — 1| and applying the conjugation transformation 
g(x) = u^ 1 oto u{x) with u(0) = 0, u{\) = 1, u'{x) > for all x G [0,1] and with the additional property 
u(x) = 1 — u(l — x) (antisymmetric with respect x ma x — \) we can span the space of all doubly symmetric maps 
(symmetric maps and symmetric invariant densities). The invariant density of the transformed map g{x) is then: 



du(x) 
dx 



To arrive in an analogous way at a nonsymmetric density while the map remains symmetric we have 



to use a generalized conjugation transformation: 

f(x) = U- 1 o t o u{x) ; U{x) =u{x)+v{x) (8) 

For symmetric f(x) the function v(x) must also be symmetric: v(x) = v(l — x). The two constraints, i.e. that 
the invariant density of the resulting map f(x) is a prescribed function p{x) and that f{x) is symmetric, determine 

X 

uniquely the transformation functions: u(x),v(x). Defining the symmetric p+(x) — \ j {p{t) — p(l — t)) dt and 

o 

X 

antisymmetric parts p-(x) = ^ J{p{t) + p(l — t)) dt (with p(x) = p+(x) + p-(x)) of the invariant measure /i(x) we 

o 

find: u(x) = p~(x) ; v(x) = p+(x). Using the fact that: p~(x) = ^(p(x) + 1 — /i(l — x)) we get the following 
expression for the desired map f(x): 

f(x) = p- 1 [l-\p(x)-p(l-x)\] (9) 

Comparing equations (||) and (^) we see that they are compatible provided that hf(x) = 1 — x, which corresponds to 
the correct choice of hf in (^) for a symmetric solution f(x). 

We can now extend the conjugation transformation approach in order to determine also nonsymmetric maps with 
a given invariant density. For the derivation of eq.(^) we have used a decomposition of the invariant measure in a 
symmetric and an antisymmetric part (with respect to x max ). This decomposition leads us to a symmetric solution. A 
general decomposition p+ (x) , /i_ (x) obeying p.+ (x) + p-{x) = p(x) is given in the following. In the case of unimodal 
maps we introduce the following components of the measure p(x): 

X 

V+(*) = l [ ( P (t) + H' f (t)p(H f (t))) dt 



o 



M- 0*0 = \ j (Pit) - H' s {t)p(H f (t))) dt ( 1.0) 



with Hf(x) given in eq.(Q). Using the generalized conjugation transformation of eq.(||) with u(x) — p-(x) and 
v(x) — p+(x) we get the solution of eq.(||) to the IFPP. Thus conjugation transformations of the form of eq. (gj) with 
u(x),v(x) chosen as in eq.(nOh span the entire space of unimodal maps with a given invariant density p(x). 



III. AN EXTENSIVE CLASS OF MAPS WITH THE BETA-DISTRIBUTION AS INVARIANT DENSITY 

The method to construct an one-dimensional unimodal map with a prescribed invariant density, as proposed in 
the previous section, applies for all possible forms of the invariant measure p(x) and the invariant density p(x), 
respectively. Nevertheless we want to consider a special class of particularly relevant invariant densities p(x), namely 
the two-parametric Beta-distribution: 

B(a + l,p+ 1) 

These densities are important due to their widespread occurence and applicability (for example in turbulent reacting 
flows, diffusion in disordered media, image classification see refs. [p0|) as well as their simple form. The special case 
of a symmetric density (a — f3) with a symmetric solution, i.e. a map obeying (f(x) = /(l — x)), which is unique, 
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has been already discussed in the literature Rj. We call this solution the doubly symmetric Beta-map (DSB). In the 
present paper we discuss the general case of an asymmetric invariant density (a ^ 0). 

We focus on two solutions of the IFPP with this general asymmetric invariant density, both of which are identical 
to the DSB in the case a — 0: 

• The first solution is symmetric (f(x) = f(l — x)). We call it the symmetric Beta-map (SB). 

• The second solution obeys the following constraint: 

f{ x h) = f{xR.) for x L , x R with p{x L ) = p{x R ), (12) 

which is only possible if a ■ > 0. We call this solution the Special Asymmetric Beta-map (SAB). In Fig.l we 
show an example of such a map for the case a = —0.9,(3 = —0.6. The corresponding invariant density is also 
presented and the notation regarding the imposed constraint of eq.(|l2|) is explained. 

In the following we consider only examples where both a and (3 possess the same sign, i.e. the above criterion is 
fulfilcd. These two types of solutions correspond of course to two different choices for the form of the function hf(x). 
For the case of the SB we have hf(x) = 1 — x. The second case is more delicate. It corresponds to an implicitly 
defined hf{x). From the condition in eq.(|T2"|) and eq. (|llj) we find: 

(1 - x R )xl = (1 - x L )xl (13) 

where 7 = %. Eq.jl^) defines hf and can be solved parametrically. Using A = ^ we get: 

^ = A l_ A 7+i ; XR = l _ A7+i (14) 

with < A < 1. Let us now apply the theoretical considerations of section II to derive the expressions for the two 
types of Beta-maps refered above. The invariant measure corresponding to the density ( Jl~T| ) is: 

ldtt a {l-tf 
( s _ _ B{a + l,(3+l,x) 

^ B{X} ~ B(a + 1,0 + 1) ~ B(a+ 1,0 + 1) [ ' 

where B(a, b, z) is the incomplete beta function Using eq.(|^) we can determine implicitly the left branch of the 
solution of the IFPP: 

B(a + 1,0+ 1, f L (x)) = B(a + l,0 + 1, x) -B(a + 1,0 + 1, h f (x)) + B(a + 1,0 + 1) (16) 

and a similar expression results also for f R (x): 

B{a + 1,0+1, f R (x)) =B(a+l,0+ 1, hj\x)) - B{a + 1,0+ 1, x) +B(a+1,0+ 1) (17) 

To get the two special solutions refered above we have to substitute in eqs.((l^,|l^) the suitable hf. As we already 
mentioned for the SB we choose hf(x) ~ h^ 1 {x) = 1 — x. Taking into account this special form of hf and the fact 
that the eqs.( ^6| , p7| ) are valid for x < x max and x > x max , respectively (for SB x max = i), we can recast eqs.(|l(], |l7| ) 
into a single equation: 

B(a + 1,0+ 1, ysB) = B(a + 1, + 1) - \B(a + 1, + 1, x) + B{0 + 1, a + 1, x) - B(a + 1,0 + 1)| xe [0, 1] (18) 

with ysB = fL,SB(x) for < x < x max and ysB — fR.,SB{x) for x ma x < x <1. To calculate the value of the SB for 
a given x we have to solve eq.([l8|) numerically using a standard root finding routine. It turns out that the bisection 
method [^2| is well suited for a fast and accurate determination of ysB- In Fig. 2a we present 4 characteristic members 
of the class of symmetric beta- maps calculated numerically using ( |l8| ) . 

Let us now return to the SAB case. The calculation of the map ysAB is a little more tedious. For a given x we use 
eq.(|l^) to find numerically a;^ (il) if a; < x max (if x > x max ). Note that for the SAB x max = ^xg- The numerical 
solution to eq.(|l3l) can be easily obtained using the bisection method. Then we apply eq.(||) to determine ysAB ■ 

r^v-li fl_i_i „ \ -/ B(a+1,0 +1, x)- B(a+ 1,0+1, x R )+B(a+ 1,0+1) ; < x < x max . . 
tS{a + L,p + i,y S AB) - I B {a + l,0+l,x L )-B(a+l,0 + l,x)+B(a+l,0+l) ; x max < x < 1 [i ^> 



4 



Again we have to find the root ysAB of eq.([L9|) numerically using the bisection method. We emphasize that the 
bisection procedure in solving eq.(]l^) could be circumvented in the case that the inverse of the incomplete beta 
function B~ 1 (a, b, z) is given explicitly. The calculated SAB maps, for the same choice of the values of a and (3 as in 
Fig. 2a, are presented in Fig. 2b. 

For both types of maps SB and SAB a dynamical route to the fully chaotic and ergodic state can be constructed 
applying the method described in 0. Using the numerical solutions to eqs. (|l^ , [l9| ) we calculated trajectories for the 
corresponding maps for different values of the parameters a and (3. Averaging over a long trajectory of each map we 
reproduced with a good agreement the invariant densities given as the input for the construction of these maps. We 
mention that, due to computational accuracy reasons, for a, (3 values close to —1 it is necessary to calculate 1 — y 
instead of y in that range of x for which y(x) ~ 1. 

Finally let us introduce the concept of the complementary beta-map, CB. Two maps y{x) and y{x) belonging to 
SB or SAB are called complementary, if their invariant densities p(x) and p(x) are transformed into each other by 
changing x to 1 — x. Considering the Frobcnius-Perron equation, which is the starting point for the construction of 
both SB and SAB maps, the following relation can be derived: 

B [a + 1, P + 1, y{x)] =B[P+l,a+l,y(l- x)} (20) 

Complementary maps agree with respect to certain dynamical features, e.g. the Lyapunov-exponent, as shown below. 



IV. ANALYTICAL AND DYNAMICAL PROPERTIES OF THE SB AND SAB 



Many features of the dynamics generated by a discrete map do not depend on the dynamical law as a whole, but 
only on its local behaviour in certain regions of the phase space. In particular for the case of unimodal maps the order 
of the maximum turns out to be responsible for some interesting quantities characterising the underlying dynamics 
like the singularities of the invariant density at x = and x = 1, the decay constant of the autocorrelation and the 
scaling parameter of the Feigenbaum bifurcation route |23[ . The order u m of a map at its maximum x max , is defined 
as the exponent in the approximation: 

f(x) im n )" m for x w z max (21) 

In an analogous way we can define the order of the map at the end points: x = 0, 1. In the following we will determine 
the corresponding exponents u m , v\ for the SB and SAB. 

(a) The SB case 

One can easily utilise the Frobenius-Perron-equation 

dy _ x a {l-xY + x l3 (l-x) a 



dx y a {l - y)P 

where y — f(x) is the SB-map. For x — > x max = i the nominator on the rhs is finite which yields — ™~ 

1 



(22) 



Assuming v m .SB to be the order of the map in the maximum, i.e. (1 — y) oc (x — |) m ' SB one gets 



dx °^ y a (l-y) 



v m ,sB = -ttt (23) 
p + 1 

Expanding eq.([l8|) around the points x — and x — 1 and taking into account that ysB — > for x — > or x — > 1 we 
find the exponents vq,sb and vi.sb to be: 

min(a,/3) + l 

vq.sb = vi.sb = — ; (24) 

a + 1 

(b) The SAB case 

Determining the order of the SAB at its maximum x = turns out to be slightly more complicated compared to the 
case of the SB due to the implicit form oihf. The derivative ^| of the SAB y — f(x) is given by the Frobenius-Perron 



equation. In terms of the parametric solution xr(X), Xl(X) of eq.(14) we get 
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dy(x) = | J£V^f , 0< x< x max 

"~ X I ~^d\~ I ~it ' x max < X < 1 



(25) 



The limit (x — x m ax) — * corresponds to the limit of the above equation for A — > 1. Higher derivatives ^-^ are 
obtained by differentiating ^ n -i with respect to A and subsequent dividing by ^§fi- or ^rfi-, respectively. The leading 
behaviour of the left or right derivative ^— %\l,r for y — > 1 (which corresponds to (a; — a; max ) — * 0) is given by: 



lm - (26) 



v ciA / 



A straightforward but tedious calculations shows that oc ( - 1 _ y ^„g + („_i) The exponent of (1 — y) vanishes for a 
value 

n = v m ,SAB = (27) 
which gives the order v m ,SAB of the SAB at the maximum. Expanding eq.(|l9|) around the endpoints x — 0, 1 we get: 

Vo,sab = mm(l, — — — -) ; i/ x S as = mm(-, — — ) (28) 
p(l + a) a 1 + a 

Acoording to the behaviour of the SB and the SAB near the characteristic points 0, l,x max we can classify the 
forms of maps which we obtain for different values of a and (3 as follows: 

• The region x ~ x max : The order of the map at x max , given by v m = , is a common quantity for all beta- 
maps. Maps with [3 < are differentiable and therefore flat at x max . With [3 approaching — 1, the maps become 
increasingly flatter. (3 > results in a cusp- like shape of the map near x rnax . 

• The zero points and 1: 

As far as the form of the maps near the points 1 and is concerned, four limiting cases with certain values for 
a and (3 can distinguished. 

- a«0, 1: 

The maps belonging to both SAB and SB show divergent derivatives /'(0) — ► oo, /'(l) — ► — oo. The 
map itself has a convex shape. As the maximum of the SAB is in [0, |], the right branch of the map, 
containing the unstable fixed point, is more shallow than that of the corresponding SB map. 

- awl, /3»0: 

At x = 0, both the SB as well as the SAB map develop a slope /'(0) — ► 1 which enables them to generate 
intermittent dynamics. Since the SB are symmetric we have /'(l) — ► 1. The derivatives of the SAB maps 
however show divergent behaviour, f'(l) — ► — oo 

The SB and SAB maps are identical to the DSB for a = (3. The derivatives of the maps /'(0) and /'(l) take 
on finite values. 

- a = 13 fa 0: The limit a = /3 — > of the DSB is the tent map /'(0) — ► 2, /'(l) — > -2. 

- a = /3« 1: 

The values of /'(0) = /'(l) are still finite, but increase unbounded in the limit a = (3 — ► 1. The map is 
extremely flat about the maximum and takes on values very close to 1 in nearly the whole interval [0, 1]. 

The dynamics generated by both the SAB and the SB result in singularities of the invariant density at x = 0,1 
for — 1 < a, (3 < 0. Although the strengths of these singularities can be directly read from eq.(|TTl) it is of interest 
to explore how they are determined through the order of the map at the maximum and at the endpoints 0, 1. The 
strength of the singularity at x — 1 is determined by the order of the maximum only. As the order of the map for 
both SB and SAB at its maximum is given by v m — -Jp^ , the expansion of the map in this region reads 

1 - y ps |1 - x max \ Vm = |1 - x max \T& (29) 
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Near x max the invariant measure is a smooth function. Using eq.(|l]) we can therefore write p(y) \dy\ oc \dx\. The 
dependence of the differential \dy\ on x is given by \dy\ « (x — j) 75 ^ 1 |da;|. If we now make the ansatz 
p(y) oc (1 — y)P for the invariant density at x rs 1, the Frobenius-Perron equation in orders of x reads 

(x-X max ) ,§ ^T (X-Xma,)^" 1 !^! CX |tfa| =*• Pj^J + ~ 1 = => P = & ( 30 ) 

We therefore conclude that the strength of the singularity of the invariant density at x = 1 for both types of maps 
(SB and SAB) is given implicitly by the order of the dynamic law at its maximum. 

The determination of the singularity at x — is more intricate. We will study in detail the case of the SB maps. 
The SAB maps can be treated in an analogous way. For the SB family the Frobenius-Perron equation can be written 
in the form: 

p(vsb) \dy S B\ = [p(x) + p(l - x)} \dx\, (31) 

r— , min(a,g) + l 

Using eq. (p4[) the SB at x » is approximately given by ysB ~ x "+ 1 . Let us make the ansatz p(ysB) oc y^ B 
for the singularity of the invariant density at y$B = 0. For x « p{x) and p{l — x) possess the singularities a and (3, 
respectively. Therefore the Frobenius-Perron equation reads 

- min(c«,/3) + l min(a,/3) + l n ■ / ~ a\ 

x « Ui x^+r !|dx| oc x min(a ' 0) \dx\ (32) 

^min(a,^) + l min(a,/3) + l . 

==> a — 1 — l = wm(a,P) 33) 

a + 1 a + 1 

Two cases have to be considered: 

• The first possibility is min(a, j3) — f3. In this case eq.(|33"|) yields a{(3 + 1) + (0 — a) — (a + l)min(a,/3) and 
therefore 

. P{a + 1) + (a - a) . 

mm(a, p) = mm a,p (34) 

a + 1 

with the solution a = a. The singularity at x = is given by the parameter a as expected. 

• The case a < (3 demands slightly more attention. The order of the map in x = and x = 1 is one, according to 
eq. ( |2^ ) . This makes the dynamics a candidate for intermittency, which is characterised by an expansion of the 
map near a marginally unstable fixed point of the form 

y(x)fax + ex z x^O (35) 

The order of the next-to leading term, z, determines essentially the statistical properties of the dynamical law 
near x = 0. In ref. |l4j| it was shown, that a symmetric map, with a maximum of order v m and of the form (|3q) 
at x = 0, generates an invariant density, which has a singularity of the order 

a = z (36) 

at x = 0. Using the implicit representation of the SB map: 

B(a + l,l3 + l,ysB) = B(a + l,p+l,x) + B(f3+l,a+l,x), x<^ (37) 
and expanding the incomplete Beta-functions around x = with the help of the formula: 

t a+l +0+2 

B(a + l,/3 + M)«— — --0— - (38) 
a -I- 1 a + 2 

we find: 
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Oj+1 x a + l X P+1 



VSB 

a+1 ~ a + 1 ^ + 1 



(39) 



It follows that: 



»'~*( 1 + J^)"~* + 7Ti*~" (40, 

Therefore the term next to the leading one is of order z = /3 — a + 1. Using eq.(|36|) and ^ m = we find 
a = /3+l — (/3 — a + 1) = a, i.e. the correct strength of the singularity at x — is reproduced. 

At the end of this section let us discuss the different kinds of dynamics which are generated by the different types 
of maps in the two classes SAB and SB. Depending on the values of a and (3, there are two eyecatching types of 
dynamical behaviour: 

• > a > [3 

The trajectories of both SB and SAB maps are characterised by an oscillatory motion around the unstable fixed 
point in the right interval of monotony of the map. The orbit leaves the region close to the fixed point, where 
the linear approximation of the map is valid, exponentially fast. This phase of the dynamics is highly correlated. 
The time which a corresponding trajectory spends in the neighbourhood of the unstable fixed point is for SAB 
maps typically larger than for SB maps. The reason for this is the fact that the derivations of the SAB maps 
are usually in a larger region closer to —1 than those of the SB maps. After leaving this regime we encounter 
an oscillatory motion of the orbit about the maximum. In this case the left and the right branch of the map are 
visited alternately. As the parts of the map visited are nearly symmetric with respect to the axis y(x) — x, a 
nearly stable orbit with period two occurs. This peculiar feature of the map is more common for the SB maps. 
The two singularities of the invariant density are build up by these oscillations. As the centre of the oscillations 
is shifted to the right of ^, the singularity at x = 1 is stronger than the one at x = 0. 

• a < (3 < 

The dynamics developed by both the SB and the SAB maps is dominated by the intermittent behaviour near 
the marginally unstable fixed point at x = 0. The trajectories show an interplay of chaotic motion and laminar 
phases in the region of x = 0. The comparison of the maps belonging to the two classes shows nicely the 
detailed balance of the system generating the same invariant density: For a SAB the exponent z in the expansion 
y(e) « x + e z , e <C 1 is much closer to 1 than for the corresponding SB map and the region of the map where this 
expansion holds is larger. Therefore the orbit tents to be captured for a longer time in this region. However, 
due to the diverging derivative of the SAB map at x = 1, the laminar region is visited less often, which partially 
cancels the above effect. 



V. LYAPUNOV EXPONENTS, AUTOCORRELATIONS AND POWER SPECTRA FOR SB AND SAB 

In this section we derive some dynamical and statistical properties of characteristic members of the SB and SAB 
classes of maps. First we consider the Lyapunov-exponents which are a measure for the degree of chaoticity. As the 
derivatives of both classes SB and SAB are obtained from the Frobenius-Perron equation, the Lyapunov-exponents 

Af^ and Af^ B , respectively can be calculated using the definition A = J p(x) In dy W 

o 



dx 



dx. The implicit form of 



y(x) prevents an analytical treatment. Therefore we have calculated Af B o and Af ^ numerically for different values 
of — 1 < a, (3 < 0, using the equations ( |l8| ) and (|l9|), respectively. Both Lyapunov-exponents show similar features as 
functions of a and (3. Most striking is the symmetry 

a SB _ K SB \SAB — \SAB f A-t \ 

i.e. the Lyapunov-exponent stays the same when switching from Beta-maps to the complementary ones, given in 
eq. (|2^) . This symmetry is even more astonishing, since it connects two very different kinds of dynamics with the same 
Lyapunov-exponent, i.e. with the same degree of chaoticity. The value of the Lyapunov exponent is maximal (ln2) 
for the doubly symmetric case a = (3 and decreases as \a — /3| increases. This can be at best seen in Fig. 3 where we 
show the Lyapunov exponent as a function of the parameters a, (3. The decrease in the degree of chaoticity for 
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the case a < (5 is caused by the frequent occurrence of almost regular intermittent intervals while in the case a > ft 
it is caused by the dominating hopping of the orbit around the unstable fixed point of the map. 

Let us now turn to the calculation of the autocorrelation function C(n) for the SB and SAB members. Again two 
types of maps in the SB and SAB classes have to be distinguished: 

• a> (3 

The correlation function is oscillating with period 2 and decays exponentially fast. This is a consequence of the 
oscillatory motion of the orbit about the unstable fixed point. The decay of the envelope is slower for maps with 
|/3 — a | large, since the corresponding map has a more shallow slope of the right interval of monotony, which 
allows long lasting oscillations. The mean decay time of the SB maps (see Fig. 4a) is much shorter compared 
to those of the SAB maps (see Fig. 4b). This is caused by the fact that the symmetry of the map restricts the 
length of the highly correlated, oscillating phases of the orbit. 

• a < (3 

The intermittent dynamics results in a correlation function, which decays algebraically in the long time regime 
according to C(n) oc n~ v . For weak intermittent dynamics the decay exponent r\ can be estimated analytically. 
For this aim we use the map: 

y(e)=e + ae z + ... 0<e<l (42) 



Assuming that the invariant density p(x) near the centre of intcrmittency is given by p(e) oc e™ < e <C 
1,-1 < a, I < z and following pj we derive the following scaling relation describing the autocorrelation function 
in the limit of distinctive intermittent behaviour, i.e. < z — 1 -C 1 and for n oo: 

C(n)wJffl + ^ " (43) 



with 



> 



and — = (z — 1) a (x + l) z 1 
n 



(44) 



where if is a constant characteristic for the system. This behaviour takes account of the long and highly 
correlated laminar sections of the orbit. The characteristic exponent rj decreases with increasing \(3 — a\ and 
thus corresponds to a slower decay of the correlation. For shorter times a transient behaviour leading into the 
algebraic decay can be seen. The correlation of systems generated by SB maps (Fig. 5a) decay faster compared 
to those SAB maps with the same values for a and (3 (Fig. 5b). The knowledge of the order of the SB maps in 
the maximum point and the zeros at x — and x = 1 allows a comparison of the characteristic exponent r) numer 
obtained numerically with the value r\ a -pprox obtained by the continuous time approximation in pTj . It shows a 
fair agreement that becomes better for maps which possess a more pronounced intermittency. In this case the 
values of a and (3 differ most. 



a 


P 


2 


Tfapprox 


numer 


-0.8 


-0.6 


1.2 


1.0 


1.09 


-0.6 


-0.4 


1.2 


2.0 


1.86 


-0.6 


-0.2 


1.4 


1.0 


0.81 


-0.8 


-0.2 


1.6 


0.33 


0.33 



TABLE I. Decay constants rj ap prox and n n u 71?. e r 



9 



The corresponding power spectrum of both the SB and SAB maps, defined by 



shows the typical power-law scaling behaviour: 



S(f) = -== I C{r)e- 2mfT dr (45) 
v 2n 



S(f) ex (46) 



for low frequencies. For higher frequencies S(f) differs from this rule since it reflects the transient regime of the 
correlation function. This behaviour of the power spectrum in the intermittent regime is illustrated in Fig. 6a, b for 
different sets of a, (3. 



VI. CONCLUSIONS 



We have presented a general solution to the inverse Frobenius-Perron problem for the case of unimodal one- 
dimensional maps. This solution allows us to cover the entire space of fully chaotic unimodal maps with given 
invariant density. It is shown that the general solution can also be obtained by a suitably defined conjugation 
transformation. Using this general representation one can find infinitely many particular solutions. To illustrate 
our method we applied it to obtain two classes of dynamical systems possessing the two parametric beta-probability 
function as the invariant density. We calculated the Lyapunov exponents, autocorrelation functions and power spectra 
for these maps and studied the properties of the obtained solutions in a twofold manner: First we observed for given 
values of the parameters of the invariant density the changes in the dynamical and statistical properties for systems 
corresponding to different partial solutions of the Frobenius-Perron equation. We then varied the parameters of 
the invariant density remaining within a certain class of partial solutions. From our studies we conclude that two 
typical properties dominate the dynamics for unimodal beta maps. The first property is the hopping of the chaotic 
trajectory around the unstable fixed point of the map which leads to an exponential decay of correlations. The 
decay constant depends on the explicit form of the special solution for given parameters of the invariant density. 
The second possibility is the intermittent scenario: the chaotic trajectory possesses extented laminar phases in the 
neighbourhood of x — 0. The correlations in this case decay with a power-law. Theoretical investigations predict that 
the characteristic exponent in this power-law decay is determined through the singularities of the invariant density 
at the center of intermittency (x = 0) as well as the derivative of the map at this point. We tested these estimations 
and found a satisfactory agreement with our numerical calculations. 

Thus our approach to the inverse Frobenius-Perron problem allows us to determine one-dimensional maps with a 
rich variety of dynamical behaviour and enables us to perform accurate calculations for physical observables even 
in the case where the ergodic limit is computationally inaccessible (intermittency). The obtained general solution 
involves the function hf and one might therefore use this degree of freedom to design dynamical laws with more 
prescribed quantities than the invariant density. One step in this direction would be to solve the inverse Frobenius- 
Perron problem with the additional constraint of a given autocorrelation function. Using the parametrisation (^3) of 
the general solution to the Frobenius-Perron equation it should be in principle possible to design a map which is 
closest to a certain given correlation behaviour. A detailed investigation of this problem is in progress p4|. 
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FIGURE CAPTIONS 



FIG. 1. The SAB map for a = —0.9, (3 = —0.6 and the corresponding invariant density. 

FIG. 2. Characteristic members of the family of maps with the beta-distribution as invariant density: (a) 4 Maps belonging 
to the SB class, (b) The corresponding 4 maps of the SAB class. 



FIG. 3. The Lyapunov exponent A aj| a for the SAB class as a function of a, f3. 

FIG. 4. The absolute value of the autocorrelation function |C ai/ g| for different values of a,/3 and in the regime a > f3 
(exponential regime): (a) for SB maps, (b) for maps belonging to the SAB class. 

FIG. 5. The absolute value of the autocorrelation function |C a ,^| for different values of a, (3 and for a < (3 (intermittent 
regime) : (a) for SB maps, (b) for maps belonging to the SAB class. 

FIG. 6. The square of the absolute value of the power spectrum S(f) as a function of the frequency / for different values of 
a,f3 and for a < (3 (intermittent regime): (a) for SB maps, (b) for maps belonging to the SAB class. 
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